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^sj Abstract. The dc response of an electrochemical thin film, such as the separator in a micro- 

battery, is analyzed by solving the Poisson-Nernst-Planck equations, subject to boundary conditions 
appropriate for an electrolytic/galvanic cell. The model system consists of a binary electrolyte be- 
tween parallel-plate electrodes, each possessing a compact Stern layer, which mediates Faradaic 
reactions with nonlinear Butler- Volmer kinetics. Analytical results are obtained by matched asymp- 

iq , totic expansions in the limit of thin double layers and compared with full numerical solutions. The 

analysis shows that (i) decreasing the system size relative to the Debye screening length decreases 
the voltage of the cell and allows currents higher than the classical diffusion-limited current; (ii) fi- 
nite reaction rates lead to the important possibility of a reaction-limited current; (iii) the Stern-layer 
capacitance is critical for allowing the cell to achieve currents above the reaction-limited current; and 
(iv) all polarographic (current-voltage) curves tend to the same limit as reaction kinetics become fast. 
Dimensional analysis, however, shows that "fast" reactions tend to become "slow" with decreasing 
system size, so the nonlinear effects of surface polarization may dominate the dc response of thin 
films. 
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Introduction. Micro-electrochemical systems pose interesting problems for ap- 
plied mathematics because traditional "macroscopic" approximations of electroneu- 
trality and thermal equilibrium fjQ , which make the classical transport equations more 
tractable [2], break down at small scales, approaching the Debye screening length. Of 

^ ■ course, the relative importance of surface phenomena also increases with miniatur- 

ization. Micro-electrochemical systems of current interest include ion channels in 
biological membranes QJ |S] and thin-film batteries |BJ [7| |UJ ^j , which could rev- 
olutionize the design of modern electronics with distributed on-chip power sources. 
In the latter context, the internal resistance of the battery is related to the nonlinear 
current-voltage characteristics of the separator, consisting of a thin-film electrolyte 
(solid, liquid, or gel) sandwiched between flat electrodes and interfacial layers where 
Faradaic electron-transfer reactions occur |11| . Under such conditions, the internal 
resistance is unlikely to be simply constant, as is usually assumed. 

Motivated by the application to thin-film batteries, here we revisit the classical 
problem of steady conduction between parallel, flat electrodes, studied by Nernst £2] 
and Brunner (1HI Hlj a century ago. As in subsequent studies of liquid ^1 Ej and 
solid |17II18| electrolytes, we do not make Nernst's assumption of bulk electroneutral- 
ity and work instead with the Poisson-Nernst-Planck (PNP) equations, allowing for 
diffuse charge in solution QJ [2J . What distinguishes our analysis from previous work 
on current-voltage relations (or "polarographic curves" ) is the use of more realistic 
nonlinear boundary conditions describing (i) Butler- Volmer reaction kinetics and (ii) 
the surface capacitance of the compact Stern layer, as in the recent paper of Bon- 
ncfont et. al |19| . Such boundary conditions, although complicating mathematical 
analysis, generally cannot be ignored in micro-electrochemical cells, where interfaces 
play a crucial role. Diffuse-charge dynamics, which can be important for high-power 
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applications, also complicates analysis [20], but in most cases it is reasonable to as- 
sume that electrochemical thin films are in steady state, due to the short distances 
for electro-diffusion. 

We stress, however, that steady state does not imply thermal equilibrium when 
the system sustains a current, driven by an applied voltage. In this paper, we focus 
on applied voltages small enough to justify the standard boundary-layer analysis of 
the PNP equations, which yields charge densities in thermal equilibrium at leading 
order in the limit of thin double layers |21| . In a companion paper |22|. we extend 
the analysis to larger voltages, at [22| and above [21] the Nernst's diffusion-limited 
current, and show how more realistic boundary conditions affect diffuse charge, far 
from thermal equilibrium. In both cases, we obtain novel formulae for polarographic 
curves by asymptotic analysis in the limit of thin double layers, which we compare with 
numerical solutions, and we focus on a variety of dimensionless physical parameters: 
the reaction-rate constants (scaled to the typical diffusive flux) and the ratios of the 
Stern length to the Debye length to the electrode separation. 

1. Mathematical Model. Let us consider uniform conduction through a di- 
lute, binary electrolyte between parallel plate electrodes separated by a distance, 
L. Our goal is to determine the steady-state response of the electrochemical cell 
to either an applied voltage, V, or an applied current, /. Specifically, we seek the 
electric potential $(X) and the concentrations C+(X) and C^(X) of cations and 
anions in the region < X < L. The Faradaic current is driven by redox reactions 
occurring at the electrodes, and we neglect any other chemical reactions, such as dis- 
sociation/recombination in the bulk solution or hydrogen production. Since we do 
not assume electroneutrality, the region of integration extends to the point where the 
continuum approximation breaks down near each electrode, roughly a few molecules 
away. In other words, our integration region includes the "diffuse part" but not the 
"compact part" of the double layer [Tl ITT1 l2"51 l2l>] . 

1.1. Transport Equations.. In the context of dilute solution theory Q]|2], the 
governing equations for this situation are the steady PNP equations: 

(1) _|( D+ ^ +M+2+ro+ j|). 

(3) -±(ug)-l«a>-**-)F. 

where F is Faraday's constant (a mole of charge), z±, /i±, and D± are the charge 
numbers, mobilities, and diffusivities of each ionic species, respectively, and e s is the 
permittivity of the solvent, all taken to be constant in the limit of infinite dilution. 
The first two equations set divergences of the ionic fluxes to zero in order to maintain 
the steady state, and the third is Poisson's equation relating the electric potential to 
the charge density. In each flux expression, the first term represents diffusion and the 
second electromigration. The Einstein relation, fi± = D±/RT, relates mobilities and 
diffusion coefficients via the absolute temperature, T, and the universal gas constant, 
R. 

Due to the potentially large electric fields in thin films, as in interfacial double 
layers, these classical approximations could break down [Tl 1251 I2I)| . For example, the 
polarization of solvent molecules in large electric fields can lower the solvent dielectric 
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permittivity by an order of magnitude, while also affecting the diffusivity and mobility. 
The solution can also become locally so concentrated or depleted of certain ions as 
to make finite sizes and/or interactions (and thus ionic activities) important. In spite 
of these concerns, however, it is reasonable to analyze the PNP equations before 
considering more complicated transport models, especially because our focus is on 
the effect of boundary conditions. 

1.2. Electrode Boundary Conditions. Although PNP equations constitute a 
well-understood and widely accepted approximation, appropriate boundary conditions 
for them are not so clear, and drastic approximations, such as constant concentration, 
potential or surface charge (or zeta potential), are usually made, largely out of math- 
ematical convenience. On the other hand, in the context of electric circuit models for 
electrochemical cells [3DJ |27| [28] , much effort has been made to describe the nonlin- 
ear response of the electrode-electrolyte interface, while describing the bulk solution 
as a simple circuit element, such as resistor. Here, we formulate general boundary 
conditions based on classical models of the double layer ^2 EH EE] , with a unified 
description of ion transport by the PNP equations. 

Our first pair of boundary conditions sets the normal anion flux to zero at each 
electrode, 

(4) D -^f(°) + M-^C_(0)||(0) = 

(5) D -^X {L) + H-*-FC-(L)^(L) = 0, 

on the assumption that anions do not specifically adsorb onto the surfaces, which holds 
for many anions at typical metal surfaces (e.g. SO4 2 , OH~, F~). The second pair 
relates the normal cation flux to the net deposition (or dissolution) flux, or reaction- 
rate density, R(C+, A$g), which in the dilute limit is assumed to depend only on 
cation concentration and potential drop, A<I>s, across the compact part of the double 
layer, originally proposed by Stern [3()j . Following the convention in electrochemistry, 
we take A$s to be the potential of the electrode surface measured relative to the 
solution. The reference potential is chosen so that the cathode, located at X = 0, is 
at zero potential and the anode, located at X = 1, is at the applied cell voltage V. 
Therefore, we have the following two boundary conditions, 

(6) D+ ^ iQ)+fi+z+F c + (0)^(0) = R(C + (0),H0)) 

(7) - D +^W f^+FC + (L)^(L) = R {C + {L)ML) ~ V)) . 

For electrodes, it is typical to assume a balance of forward (deposition) and backward 
(dissolution) reaction rates biased by the Stern voltage with an Arrhenius temperature 
dependence, 

(8) R(C + , A* s ) = K C C + exp ( ^g^ ) - K a C M exp (^^ 

where Cm is the (constant) density of electrode metal and K c and K a are rate con- 
stants for the cathodic and anodic reactions [y. (Expressing the reaction rate in 
terms of the surface overpotential, r/s = A$s — A$g 9 , where A$^ 9 is the Stern- 
layer voltage in the absence of current, R = 0, yields the more common form of the 
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Butler- Volmer equation ^|2S|-) The Stern layer voltage contributes — zFA$$ to 
the activation energy barriers multiplied by transfer coefficients a c and a a for the 
cathodic and anodic reactions, respectively, where a c « a a fa h, for single electron 
transfer reactions [Tl 1261 1^9"] . 

Following Frumkin |31j . we apply Eq. JSJ just outside the Stern layer, in contrast 
to "macroscopic" models which postulate the Butler- Volmer equation as a purely 
empirical description of reactions between the electrode surface and the electrically 
neutral bulk solution. Physically, the Frumkin approach makes more sense since the 
activation energy barrier described by the Butler- Volmer equation actually exists at 
the atomic scale in the Stern layer, not across the entire "interface" including diffuse 
charge in solution. We are not aware of any prior analysis with the full, nonlinear 
Butler- Volmer equation as a boundary condition on the PNP equations other than 
that of Bonnefont et al. J5j . Earlier analyses by Jaffe et. al. |15l l*H)] and Itskovich 
et. a/.|17| also include electrode reactions, but only for small perturbations around 
equilibrium. 

The final pair of boundary conditions determines the electric potential by spec- 
ifying the voltage drop across the Stern layer in terms of the local electric field and 
concentrations , 

(9) - $(0) = A$ s (JL(0), C + (0), CL(0) 

(10) V-*{L) = M 3 l~(L),C+(L),C-(L) 

In macroscopic electrochemistry, these boundary conditions are usually replaced by 
the assuption of electroneutrality (which eliminates the need to solve Poisson's equa- 
tion) or by simple Dirichlet boundary conditions on the potential [lj. In colloidal 
science, they are likewise replaced by simple boundary conditions of constant surface 
charge (or zeta potential) [221 EH]- Here, we incorporate more realistic properties of 
the interface as follows: Neglecting the specific adsorption of anions, the Stern layer 
acts as a nonlinear capacitor in series with the diffuse layer. Grahame's celebrated 
electrocapillary measurements jSU EH] suggest that (i) the Stern layer capacitance, 
Cs, is roughly independent of concentration, depending mainly on the (variable) total 
charge, a, 

(11) < A ^- ' 



da Cs(a) 1 

and (ii) dilute solution theory accurately describes the capacitance, Cd (c, C+ ) , of 
the diffuse layer, at least when the charge and current are small enough to be well- 
described by Poisson-Boltzmann theory (as derived below). Using Gauss' law, the 
surface charge density can be expressed in terms of the normal electric field, a = 
—esd^/dX, where es is an effective permittivity of the compact layer. Therefore, 
integrating Eq. i|ll|l . Grahame's model corresponds to the assumption, 

-tsd^/dx . 

(12) A$ s = ' 



Cs(°Y 



which determines how the voltage across the compact layer (relative to the point of 
zero charge for which A<i>s = 0) varies as the two capacitors become charged. The 
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function, Cs(cr), should be fit to experimental or theoretical electrocapillary curves 
at large concentrations (since l/C tota i — 1/Gx> + 1/Cs ~ 1/Cs m that case). 

The simplest model that captures this interplay between the compact and diffuse 
layers is the Stern model ^2 EH- which assumes the capacitance of the compact 
layer, Cs, to be constant [30. While more complicated models for the compact 
layer have been proposed [SEl EH E3 the Stern model suffices for our purposes, 
because it allows us to describe surface capacitance easily in the context of our model 
of Faradaic reactions. Following Itskovich et al. ^7] and Bonncfont et al. ^!5]j let 
us introduce an effective width, As, for the compact layer, As = es/Cs, so that 
Eq. (|12|l reduces to a linear extrapolation of the potential across the compact layer, 
— A$s = Xsd^/dx. Substituting this expression into Eqs. © and I|1U|) yields two 
Robin boundary conditions, 

(13) $ (0 )-As||(0) = 

(14) $(L) + \s^(L) = V, 

completing a set of six boundary conditions for our three second-order differential 
equations. Physically, the Stern layer, as an effective solvation shell for the electrode, 
is only a few molecules wide, so it is best to think of As as simply a measure of the 
capacitance of the Stern layer. More generally, the same boundary condition could 
also describe a thin dielectric layer on the electrode |381 1391 130) . e.g. arising from 
surface contamination or a passivating monolayer. 

Note that since the anion flux is zero, the current passing through the cell is 
proportional to the cation flux (everywhere in the cell, since it is constant), 

. / dC+ d<$>\ 

(15) I = Z+ FA[D + ^± + » +Z+ FC + -), 

where A is the electrode area and a current flow towards the cathode (x — 0) is 
taken to be positive. Under potentiostatic conditions, the cell voltage V is given, 
and the steady-state polarization curve I{V) is determined by solving the equations. 
Conversely, under galvanostatic conditions, / is fixed, and we solve for V(I). 

1.3. An Integral Constraint. As formulated above, the boundary value prob- 
lem is not well-posed. Since the anion flux is constant throughout the cell according to 
Eq. @ , the two anion flux boundary conditions are degenerate, leaving one constant 
of integration undetermined. This is not surprising as we have omitted one crucial 
physical parameter, the total number of anions. More precisely, because anions do 
not react, we must specify their total number, which remains constant as the steady 
state is reached. This corresponds to the constraint, 

(16) ~ / C_(X)dX = C ref) 

L Jo 

where C re f is the initial concentration of anions. 

Note that the total number of cations (and hence the total charge) is not known 
a priori because the removal of cations at the cathode and the injection of cations at 
the anode may significantly alter their total number. This may seem counter-intuitive 
since we are accustomed to assuming that we know the total cation concentration 
at all times based on the original molarity of the solution, but this "macroscopic" 
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thinking does not apply when the physics at the microscopic level are explicitly being 
studied (e.g. diffuse charge layers or micro-electrochemical systems). Mathematically, 
the reaction boundary conditions at the electrodes Eqs. © and J7J) are sufficient to 
determine the total cation concentration (and total charge), as long as the total anion 
concentration is specified. 

1.4. Dimensionless Formulation. To facilitate our analysis, we formulate the 
problem in dimensionless form. For simplicity we also assume that the electrolyte is 
symmetric, z + = — z_ = z, which does not qualitatively affect any of our conclusions 
as long as z + /|z_| is not too different from 1 (which holds for most simple, aqueous 
electrolytes). Scaling the basic variables as follows, 



(17) x = X/L, 


c±(x) = C±(xL)/C ref 


Eqs. (d-© become 




(18) 


d 2 c + d ( d(j) 
dx 2 dx \ dx 


(19) 


d 2 c^ d ( dej) 
dx 2 dx \ dx 


(20) 


2 d 2 4> l 

~ € d^ = 2 (C+ 



<j>(x) 



zF${xL) 
KT ' 



dc+ 

dx 


d(/) 
dx 


= 4i 


dc_ 

dx 


dx 


= 0, 



where e = Xd/L is the ratio of the Debye screening length \e> = J 2z^f?c — ^° ^ e 
distance between electrodes. The Debye length is typically on the order of nanometers, 
so e is always extremely small for macroscopic electrochemical cells. This situation 
changes, however, as L or C re f is decreased, and in the case of nano-electrochemical 
systems e could be as large as 10. 

The two flux equations are easily integrated, using Eq. Q to evaluate one constant 
and leaving the other constant expressed in terms of the current via Eq. (|15fl > 

(21) 

(22) 

where we have defined a dimensionless current density, i = I /Id, scaled to the Nernst's 
diffusion-limited current density (see Sect ion |2~T|) . Id = AzFD + C re fA/L. Since diffuse 
charge is of primary interest here, it is convenient to introduce 

(23) c=-(c + + c-) and p=-(c+-C-), 

the average concentration of ions and (half) the charge density, respectively, which 
leaves us with a coupled set of one second-order and two first-order differential equa- 
tions, 

fr, A \ dc d(f> 

dx dx 
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Nondimensionalizing the boundary conditions (remembering that there is an extra 
condition needed to determine the current-voltage relation, i(v) or v(i)), we obtain: 

(27) 0(0) - fc^O) = 0, 

ax 

(28) 0(l)+<jAl)=«, 

ax 

(29) kc[c(0) + p(0)}e a ^ (0) - i r e- a ^ (0) = i, 

(30) -Jfcc[c(l) + p{l)]e a " Wl) - v) + i r er a ^ {1) - v) = i, 

(31) f [c(x) - p{x)]dx = 1, 

Jo 

where 

32 fc„ = , i r = — — , v = , and o = . 

V ' c 4D+ ' 4D+C ref ' iff ' X D 

Keep in mind that the dimcnsionless rate constants decrease with system size, so that 
"fast reactions" (k c ,i r ^> 1) may become "slow reactions" (k c ,i r — 0(1)) as L is 
reduced to the micron or submicron scale. 

It is important to note that we have scaled the effective Stern layer width, As, 
with the Debye screening length, Ad, rather than the electrode separation L, thus 
introducing the factor e = Xd/L in Eqs. (|27|l and (|28|l . This choice is important for 
our asymptotic analysis of the limit e — > at fixed 8, which is intended to describe 
situations in which L is much larger than both Xs and Ad. Without it, our analysis 
would assume that as e — ► the Stern layer becomes infinitely wide compared to 
the diffuse layer, even though it is mainly the macroscopic electrode separation which 
varies. The limit of very small Stern layer capacitance, which amounts to the Helmoltz 
model of the double layer 0J, is best studied by letting 8 — > oo after e — > 0. In 
contrast, because e and 8 would both be small, the limit of very large Stern layer 
capacitance can be studied by simply letting 8 = 0, yielding the Dirichlct boundary 
conditions, 

(33) 0(0) = , 0(1) = 1, 

of the Gouy-Chapman model of the double layer |421 143) . In our work, we shall 
consider both limits, starting with the assumption that 8 = 0(1), which corresponds 
to the Gouy-Chapman-Stern model of the double layer [TT1I21)] . 

For one-dimensional problems, galvanostatic conditions are more mathematically 
convenient than potentiostatic conditions. In the former case, i is given, and v(i) is 
easily obtained from the Stern boundary condition at the anode Eq. 1)28)1 . In the 
latter case, however, v is specified and i(v) must be determined self-consistently to 
satisfy Eq. 1)28)) . Therefore, we shall assume that the current i is specified and solve 
Eqs. 1)24)1 - 1)26)1 subject to the boundary conditions Eqs. ()29)l and 1)30)) and the integral 
constraint Eq. 1)31)1 . 

2. Boundary-layer Analysis. In this section, we briefly review the classi- 
cal asymptotic analysis of the PNP equations, pioneered independently by Cher- 
nenko 03j, Newman [21], and MacGillivray (3Hji which involves boundary layers of 
width, e (corresponding to diffuse-charge layers of dimensional width, Ad). As dis- 
cussed below, the classical asymptotics breaks down at large currents approaching 



8 MARTIN Z. BAZANT, KEVIN T. CHU, AND B. J. BAYLY 

diffusion limitation. Unlike most previous authors, who assume either a fixed po- 
tential 03 EH] or fixed interfacial charge |2J a t an isolated electrode or fixed con- 
centrations at cell boundaries with ion-permeable membranes 0] 02] , we solve for 
the response of a complete, two-electrode galvanic cell with boundary conditions for 
Faradaic reactions and Stern-layer capacitance. 

Throughout this section, the reader may refer to Figure ^ which compares the 
uniform asymptotic solutions derived below to numerical solutions at several values of 
e. These figures illustrate the structure of the field variables in the cell as well as give 
an indication of the quality of the asymptotic solutions. The numerical solutions are 
obtained by a straightforward iterative spectral method, described in a companion 
paper [22] , 

2.1. Electroneutrality in the Bulk Solution . The most fundamental ap- 
proximation in electrochemistry is that of bulk electroneutrality \f\. As first em- 
phasized by Newman J5J, however, this does not mean that the charge density is 
vanishing or unimportant, but rather that over macroscopic distances the charge den- 
sity is small compared to the total concentration, \C+ — C_| <C C+ + C_, or, in our 
dimensionless notation, \p\ <C c. Mathematically, the "macroscopic limit" corresponds 
to the limit e = Xd/L — ► 0. The electroneutral solution is just the leading order so- 
lution when asymptotic series of the form f(x) = f^°'(x) + ef^'(x) + e 2 f^ 2 \x) + . . ., 
are substituted for the field variables in Eqs. l|24 fl -l|26 fl . 

Carrying out these substitutions and collecting terms with like powers of e, we 
obtain a hierarchy of differential equations for the expansion functions. At 0(1) we 
have, 

(34) ^ = 2: , -<#»E<n = 2i , pW = 

ax 

where the bar accent indicates that these expansions are valid in the "bulk region" 
e <C x <C 1 — e (or Ac cXcL-Ad). Integating these equations, we obtain the 
leading order bulk solution: 

(35) c (0) (x) = c + 2ia; 

(36) P(s) 



-1 



x + c D /2i 
(37) ^(o) (x)=0o + log (\ + ^ 

where the constants of integration c and <j) are the values of the bulk concentration 
and potential extrapolated to the cathode surface at x = 0. Note that, despite quasi- 
electroneutrality, the electrostatic potential does not satisfy Laplace's equation at 
leading order in the bulk, as emphasized by Levich 0Sj and Newman pp. Noting the 
presence of e 2 in Eq. Ij26|l for the dimensionless potential, it is clear that a negligible 
charge density, p = 0(e 2 ), is perfectly consistent with a nonvanishing Laplacian of 
the potential. More precisely, we have, 

(38 , *»w-*a,) ' 



dx 2 v ; (x + c /2i) 2 



at 0(e 2 ) in Eq. ^gj. 
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Fig. 1. Numerical solutions (solid lines) compared with the leading order uniformly-valid 
approximations (dashed lines) given by Eg. 15H6 - 161V for the dimensionless potential, </>(x), electric 
field, E(x), concentration, c(x), and charge density, p(x) for the case j = 0.9, k c = 10, i r = 10, 
5 = 0, and e = 0.001,0.01,0.1. Linear scales on the left show the entire cell, while log scales on 
the left zoom in on the cathodic region. Note that in the concentration and charge density profiles, 
the numerical and asymptotic solutions are barely distinguishable for e < 0.01. 

The integral constaint, Eq. (|31|l . can be used to evaluate the constant c . If we 
assume that c_ in the boundary layers does not diverge as e — > 0, then they only 
contribute O(e) to the total anion number. Therefore, 



(39) 



c-(x)dx= I c iQ) (x)dx + O(e) = c + i + O(e), 



which implies that c = 1 — l. 
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At leading order in the bulk, we have recovered the classical theory dating back 
a century to Nernst ^3 1131 HI] . The solution is electrically neutral with a linear 
concentration profile whose slope is proportional to the current. This approximation 
leads to one of the fundamental concepts in electrochemistry, that there exists a 
"limiting current" , l = 1 , or 

, s 4zFD+C ref A 
(40) I = h = jf**—> 

corresponding to zero concentration at the cathode, c = 1 — l = 0. The current is 
limited by the maximum rate of mass transfer allowed by diffusion, and larger currents 
would lead to unphysical and mathematically inconsistent negative concentrations (see 
Appendix ^J . 

Examination of the electric field exposes the same limitation on the current: a 
singularity exists at i = 1 that blocks larger currents from being attained. The leading 
order bulk approximation to the electric field, E^ — l/[x + (1 — i)/2i], diverges near 
the cathode like 1/x in the limit l — ► 1. This would imply that the cell voltage 
V (calculated below) diverges as i — > 1, thus providing a satisfactory theory of the 
limiting current since an infinite voltage would be necessary to exceed (or even attain) 
it. Unfortunately, this classical picture due to Nernst ^21; based on passing to the 
singular limit e = 0, is not valid for any finite value of e because the solution is, in 
general, unable to satisfy all of the boundary conditions. 

2.2. Diffuse Charge Layers in Thermal Equilibrium. We now derive the 
leading order description of the boundary layers in the standard way |211 151 051 IT]. 

using Eqs. (|24Jl -125 |l . The singular perturbation in Eq. 1261) can be eliminated with 
the rescaling y — x/e indicating that the boundary layer at x = has a width O(e). 
In terms of this inner variable, the governing equations in the cathode boundary layer 
are 

(41) 

(42) 

(43) 

where e now appears as a regular perturbation since solutions satisfying the cathode 
boundary conditions and the matching conditions still exist when e = 0. At the anode, 
the appropriate inner variable is y — (1 — x)/e, and the equations are the same as 
above except that i is replaced with — l, since current is leaving the anode layer, while 
it entering the cathode boundary layer. 

Expanding the cathode boundary layer fields (indicated by the check accent) as 
asymptotic series in powers of e, we obtain at leading order, 

(44) ^ + /»>f!U 

dy dy 

(45) ^ + ^ =0 . 

dy dy 

Using Eq. I|23|) to rewrite these equations in terms of c+ and c_ , we find that the flux 



dc dcj) 
dy dy 


= 2ie 


dp dcj) 
dy dy 


= 2ie 


d 2 c)> _ 
dy 2 


= P, 
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of each ionic species in the boundary layer is zero at leading order: 

(A^ d4 0) , _(0)4 (0) n 

(46) ^T ±C± ^r =0 " 

While this equation appears to contradict the fact that the current is nonzero, the 
paradox is resolved in the same way as electroneutrality is reconciled with a non- 
harmonic potential in the bulk region: tiny fluctuations about the boundary layer 
equilibrium concentration profiles at O(e) are amplified by a scaling factor of 1/e to 
sustain the O(l) current. Thus, the leading order contribution to the current in the 
boundary layer is 

Uc ( l ] s(1) #(°) M 0*>\ _ dc^ _ (1) 4 (0) ^ .(O).-' 11 ' 



\ dx dx dx J dy dy dy 

Integrating Eq. (|46|l and matching with the bulk, we find that the leading order ionic 
concentrations are Boltzmann equilibrium distributions 1 : 

(48) i )(y)=c o e ± ^-^ ( ° ) ^ 

where c = 1 — j and cf> = ft ' (0) are obtained by matching with the solution in the 
bulk. Note that the Boltzmann distribution arises not from an assumption of thermal 
equilibrium in the boundary layer but as the leading order concentration distribution, 
even in the presence of a non-negligible 0(1) current. 

The general leading-order solution was first derived by Gouy 021 an d Chap- 
man 03] and appears in numerous books Ql EM G21 EH an d recent papers |5HT5): 

(49) c(%)=c cosr# -<£(%)] 

(50) p(°>(y)= Co sinl# -^%)] 

(51) __ =2 ^ smn ^ _ j 

(52) ft°\y) = O + 4tanh- 1 ( 7o e-^) , 

where 7 = tanh(Co/4) and Co = < -°- ) (O) — cf> is the leading order "zeta potential" 
across the cathodic diffuse layer, which plays a central role in electrokinetic phenom- 
ena (221 US]- Note that the magnitude of the diffuse layer electric field scales as 1/e 
as illustrated in Figure ^ 

The value of ft°^ (0) hidden in the zeta potential Co is determined by the Stern 
boundary condition, Eq. i|27|) . If 5 — (Gouy-Chapman model), then ft°^(0) = 0, or 
Co = — 4>o, which means that the entire voltage drop 4> across the cathodic double 
layer occurs in the diffuse layer. If 5 = oo (Helmholtz model), then <^°'(0) = 4> , or 
Co = 0, in which case the Stern layer carries all the double layer voltage. For finite 
S > (Stern model), Co is obtained in terms of <f> by solving a trancendental algebraic 
equation, 

(53) - Co = 2<5V5o"sinh(C/2) + 4> , 



1 The expresssion for energy in the Boltzmann equilibrium distribution includes only the energy 
due to electrostatic interactions. "Chemical" contributions to the energy are neglected. 
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which can be linearized about the two limiting cases and solved for £ , 
b - 26y/c;smh(<p /2) if 6 < <p /2^/c;smh(fa/2) 



(54) C ° ' " 1 fa/5^T if 6 » o /2^ 



Note that if <^> <SC 1 , then — £ rs rzf~/ — ^ s a reasonable approximation for any value of 
<5 > 0. Finally, we solve for O by applying the Butler- Volmer rate equation, Eq. (|25|) . 
which yields a transcendental algebraic equation for <f> : 

(55) fc cCoe -^+«^+« _ lre -M< o +0o) = L 

Simultaneously solving the pair of equations (|53[) and l|55|) exactly is not possible in 
general, but below we will analyze various limiting cases. 

In the anodic boundary layer, we find the same set of equations as Eqs. (|4T]l - fT3l) 
except that i is replaced by — i. Therefore, since the fields do not depend on i at 
leading order, the anodic boundary layer has the same structure but with different 
constants of integration. Thus, we find that the leading order description of the anodic 
boundary layer is given by Eqs (|49|I - H52|) with c , <p , j , and £ replaced by different 
constants c\, fa, 71, and £1, respectively. Moreover, it is straightforward to show that 
Cl = g<°) (1) - 1 + 1 and fa = fa + log (i±i) . 

The leading-order anodic zeta potential, £1, and potential drop across the entire 
anodic double layer, v — fa, are found by solving another pair of trancendental alge- 
braic equation resulting from the anode Stern and Butler- Volmer boundary conditions, 
Eqs. (Pl> and pi|). 

(56) - Ci = 25 V^T sinh(Ci /2) + fa - v 

(57) 1 = -fc c cie- Cl+Q ^ Cl+ * 1 -^ + i re -«.Ki+*-«). 

As before, the Gouy-Chapman and Helmholtz limits are £i = v — fa and £1 = 0, 
respectively, and for small voltages (or currents) the approximation £i rs (v — fa)/(l + 
Sy/ci) is valid for all 5 > 0. 

2.3. Leading Order Uniformly-valid Approximations. We obtain asymp- 
totic approximations that are uniformly valid across the cell by adding the bulk and 
boundary layer approximations and subtracting the overlapping parts: 

(58) c(x) = $°\x/e) - c ] + c^(x) + [c<°>((1 - x)/e) - Cl ] + 0(e), 

(59) p(x) = p (0) (x/e) + e 2 p^(x) + p^((l - x)/e) + 0(e), 

60 E(x) = — ^— (x/e) + -^— x ^- 1 - 1 /e + 0(e), 

e ay ax e ay 

(61) fax) = [ft°\x/e) - fa] + 0<°>(aO + [0 (o) ((l - *)/e) - fc] + 0(e). 

Note that we have kept the 0(e 2 ) term in the charge density since it is the leading 
order contribution in the bulk region and is easily computed from Eq. (|38[) . As shown 
in Figure^ the leading-order uniformly valid solutions are very accurate for e < 0.01 
(or L > IOOAd) and reasonably good for e = 0.1. Since higher-order terms are not 
analytically tractable, it seems numerical solutions must suffice for nanolayers, where 
e rs 1, or else other limits of various parameters must be considered, as below. 

The discrepancy in electric potential profile at large e in Figure^is particularly in- 
teresting because it arises from a constraint on the total potential drop across the cell. 
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To understand the origin of this voltage constraint, recall that the total cell voltage 
is determined by the current density flowing through the cell (via the voltage-current 
relationship). While e is technically a parameter in the voltage-current relationship, 
a leading-order analysis does not capture the e dependence. Thus, the leading-order 
cell voltage must be the same for all e which is what we observe in figure ^ A close 
examination of the potential and electric field profiles reveals that most of the error in 
the asymptotic solution for the potential comes from an over prediction of the electric 
field strength (and therefore the potential drop) in cathode region. 

3. Polarographic Curves for Thin Double Layers, e — ► 0. The relationship 
between current and cell voltage is of primary importance in the study of any elec- 
trochemical system, so we now use the results from the previous section to calculate 
theoretical polarographic curves in several physically relevant regimes. We focus on 
the effects of the Stern capacitance and the reaction rate constants through the di- 
mensionless parameters, 6, k c and i r , with a c — a a = 1/2. For a fixed voltage, the 
mathematical results are valid in the asymptotic limit of thin double layers, e — > 0. 

3.1. Exact Results at Leading Order. Using the uniformly valid approxi- 
mation Eq. 161|) . we can write the leading order approximation for the cell voltage 
as 

(62) u = o + 2tanlT 1 (i) + (v-0 :L ). 

We can interpret this expression as a decomposition of the cell voltage into the poten- 
tial drop across the cathode, bulk, and anode layers respectively. Note the divergence 
in the bulk contribution to the cell voltage as i — ► 1, which we expect from our earlier 
analysis. In the next section, we explore analytic solutions for several limiting cases 
and compare them to exact solutions given by Eq. I|62|) with the leading-order cath- 
ode and anode diffuse layer potential drops determined implicitly by Eqns. H53|) . Q55|l. 
(|5t>p. To make plots in our figures, we use Newton iteration to solve for <f> and v — 4>i 
in this algebraic system. 

3.2. Cell Resistance at Low Current . Given the common practice of using 
linear circuit models to describe electrochemical systems |27ll28ll^T)j . it is important to 
consider the low-current regime, where the cell acts as a simple resistance, R = V/I. 
First, we compute the potential drop across the double layers. Since the procedure 
is almost identical for the two boundary layers, we focus on the calculation for the 
cathode. By writing the boundary conditions Eq. (|29|l in the standard Butler- Volmer 
form involving the exchange current and surface overpotential |26l ^ : 



(63) i = f (e 



~ a cV s 



where i c = (fc c c e~^°) a i" c and if s — A0s — A<fi e g are the cathode exchange current 
and surface overpotential, respectively. Note that the exchange current contains the 
Frumkin correction through the factor e~^° |26| . For low current densities, we expect 
the surface overpotential to be small, so we can linearize this equation to obtain 

(64) i ~ -i% 

where we have used the fact that a c + a a = 1 . Rewriting this equation in terms of 
(f)(0), we find that 

(65) 0(O)~^+<MO), 

1 
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where e q(O) is the value of 0(0) calculated from the cathode Butler- Volmer rate 
equation when there is no current flowing through the electrode. The zeta potential 
Co in the formula for the exchange current is determined by combining Eq. I|65|) with 
(I55J1 to obtain a single equation for £ : 



(66) -2<ty£sinli(C/2)~- _\ - +log ( £- ) + Co. 

Finally, to compute the total double layer potential drop we add the potential drop 
across the diffuse layer to 0(0): 

(67) O = 0(0) - Co ~ - c + log (^- 

1 \ KcCo 

A similar calculation at the anode results in 

(68) ^-^l+iogfe 

where i" = (fccCie - ^ 1 ) ° i" c and Ci is determined by the anode equivalent of Eq. (|66J1 . 
Combining these results with the potential drop across the bulk solution, we find 
that the total cell voltage is given by 

v(i) ~ 4tanh~ 1 (i) H 1 

y y y ' -,c -.a 

(69) mU+L + L 



This result gives the dimensionless resistance, r, of the electrochemical thin film as a 
function of the physical properties of the electrodes and the electrolyte. Note that the 
Stern-layer capacitance is accounted for implicitly via the calculation of the electrode 
zeta potentials. 

3.3. Simple Analytical Formulae. The exact leading-order current-voltage 
relation simplifies considerably in a variety of physically relevant limits. These ap- 
proximate formulae provide insight into the basic physics and may be useful in inter- 
pretting experimental data. 

3.3.1. The Gouy-Chapman Limit (S — ► 0). In this limit, the capacitance of 
the diffuse layer of the charged double layer is negligible compared to the capacitance 
of the compact layer. As a result, the voltage drop across the diffuse layer accounts 
for the entire potential drop across the charged double layer. Physically, this limit 
corresponds to the limits of low ionic concentration or zero ionic volume |26j . Since 
Co = — O and Ci = v — 4>i when 6 = 0, the Butler- Volmer rate equations, Eqs. I|55|l 
and 157J1 . reduce to 

(70) fc c (l - i)e~ c ° - v = i and - fc c (l + i)e~ Cl + i r = i. 
Solving for Co and Ci > we find that 

(71) Co = ln(l-i)-ln(^-^ and Ci=ln(l + i) + h/ '''' 
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GC Limit (S = 0) 

"0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

J 

Fig. 2. Exact polarographic curves (dashed lines) for varying 8 values compared to polaro- 
graphic curves for the Gouy-Chapman (8 = 0) and Helmholtz (8 = ooj limits (solid lines). Top: a 
reaction-limited cell (\ r < I) with physical parameters, k c = 0.03, i r = 0.7. Notice that above the 
reaction-limited current density, \ r , the highest cell voltages occur for 8 values near 0. Bottom: 
a diffusion-limited cell (i r > 1) with physical parameters: k c = 0.05, i,- = 1.5. In both cases, 8 
increases as the curves move upwards. 



which can be substituted into Eq. I|tj2|) to obtain 

(72) v(i) = 4tanh~ 1 (i) + 2tanh _1 (i/i r ). 

Notice that the boundary layers make a nontrivial contribution to the leading order 
cell voltage. The 2tanh~ 1 (i/i. r ) term is especially interesting because it indicates the 
existence of a reaction-limited current when i r < 1. In hindsight, it is obvious that re- 
action limited currents exist in the Gouy-Chapman limit because the reaction kinetics 
at the anode do not permit a current greater than i r . We emphasize, however, that 
the Gouy-Chapman limit is singular because there is no problem achieving current 
densities above i r for any S > (see Figure 0). For any finite S > 0, the shift of the 
anode double-layer potential drop to the Stern layer helps the dissolution reaction 
while suppressing the deposition reaction which permits the current density to rise 
greater than i r . 

Note that the cathodic and anodic boundary layers do not evenly contribute to the 
cell voltage near the limiting currents. In a diffusion-limited cell, the cathodic layer 
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makes the greater contribution because as 1 — ► 1, £ diverges while £i approaches a 
finite limit. We expect this behavior because as 1 — ► 1, the electric field only diverges 
at x = 0. However, when the cell is reaction-limited, the division of cell voltage 
between the boundary layers is reversed as i approaches the limiting current i r . Even 
the voltage drop in the bulk becomes negligible compared to Ci in the reaction-limited 
case. In this situation, the cell voltage diverges as i — > i r because the only way to 
achieve a current near i r is to drastically reduce the deposition reaction at the anode. 
In other words, the cation concentration at the anode must be made extremely small 
which requires a huge anodic zeta potential. 

3.3.2. The Helmholtz Limit (<5 — ► oo). This is the reverse of the Gouy- 
Chapman limit. Here, the capacitance of the compact layer is negligible, so the 
potential drop across the double layer resides completely in the compact layer. The 
Helmholtz limit holds for concentrated solutions or solvents with low dielectric con- 
stants and other situations where the Debye screening length becomes negligible J2E] ■ 
It also describes a thick dielectric or insulating layer on an electrode [381 I2U| . 

In the Helmholtz limit, £ = = £i, so the Butler- Volmer rate equations take the 
form 

(73) k c {l ~ i)e a ^° ~ ire-^^ =i 

(74) -k c (l + i) e a *Wi-*) + i re - Q " (01 - l,) = l. 

Solving these equations for <f> and v — </>i under the assumption of a symmetric 
electron-transfer reaction (i.e. a c = 1/2 = a a ) and substituting into the formula for 
the cell voltage, we find that 



(75) v{\) = 6tanh~ 1 (i) + 21n 



i+ y / i 2 +4i r fc c (l-i) 
-i+ v /i 2 + 4i,.fc c (l + i) 



While this expression appears to be more complicated than the one obtained for the 
Gouy-Chapman model, it is not very different when i r > 1 as can be seen in Figures|21 
and |31 In fact, the wide spread in the polarographic curves observed in Figure [5] 
requires that k c <C i r ] otherwise, all of the curves would be difficult to distinguish. 
Moreover, as we shall see in the next section, in limit of fast reactions, both models 
lead to the same expression for the cell voltage for i r > 1. On the other hand, when 
i r < 1, the two models are qualitatively very different. While the Gouy-Chapman 
model gives rise to a reaction-limited current, the Helmholtz model does not. 

3.3.3. The Fast-Reaction Limit, (i r > l,(i r ) a "(k c ) ac > i/(l -i) Q °) • The 
polarographic curves for all 8 values collapse onto each other in the limit of fast 
reaction kinetics (Figure|3J). Even the assumption of symmetry in the electron-transfer 
reaction is not required. When reaction rates are much larger than the current, the 
two reaction-rate terms in the Butler- Volmer equations, Eqs. <|55[1 and JS7J, must 
balance each other at leading order: 

(76) k c (l - !) e -Co+«c(Co+0o) _ lre - aa (Co+M M 

(77) -k c (l + 1 ) e -Ci+«c(Ci+0i-^) + lre -MO+0i-«) ~ . 

Since a c + a a = 1 for theoretical models of single electron-transfer reactions [Tl lllll2T)) . 
we can solve explicitly for (f> and v — 4>i to obtain 

(78) *-='"U<f^)J - »-*='»' M1+,) 
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Thus, for fast reaction kinetics, the leading order cell voltage is given by 

(79) v(i) = 4tanh~ 1 (i). 

Notice that this is exactly the fast reaction limit of v(i) that we find in both the 
Gouy-Chapman and Helmholtz limits. It is straightforward to check the validity of 
the assumptions made in Eqs. I|7tj|) an d (|77|l by substituting these results into the 
expresssions for the reaction rates and observing that the zeta potentials satisfy the 
bounds Co < and Ci < hi ( — _^ ) which follow from the monotonicity of £ and £i 
as functions of S. 



Polarographic Curves in Gouy-Chapman Limit 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
J 

(a) 
Polarographic Curves in Helmholtz Limit 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

J 

(b) 

Fig. 3. Polarographic curves in the (a) Gouy-Chapman limit and (b) Helmholtz limit as the 
reaction rate constants are increased (dashed lines). For these plots, the reaction rate constants 
(\ r = 1, 2, 5, andlO) increase as the curves shift towards the lower right and are related by k c = 
i r /2. It should be noted that the fast reaction limit is reached very guickly; in both plots, the curve 
closest to the fast reaction curve has a i r value of only 10. 
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4. Thick Double Layers, e = 0(1). Up to this point, we have only examined 
the current-voltage characteristics in the singular limit e — ► 0, where the current 
density cannot exceed its diffusion- limited value, i = 1. The situation changes changes 
for any finite e > 0. 

4.1. What Limiting Current?. As is clearly evident in Figure^] the cell has 
no problem breaking through the classical limiting current for e > 0. Figure 0] also 
shows that the e dependence of the polarographic curves only becomes significant at 
currents approaching the diffusion- limited current; below i « 0.5, the curves are nearly 
indistinguishable. Moreover, as e increases, upper end of the polarographic curves 
flatten out and shift downwards. This decrease in the cell voltage for large e values 
arises because the diffuse charge layers overlap and are able to interact with each other. 
More precisely, the cell has become so small (relative to the Debye screening length) 
that the electric fields from the two diffuse layers partially cancel each other out 
throughout the cell resulting in a lower total cell voltage. It should be emphasized that 
this effect is only observable because we are studying a two electrode system. Single 
electrode systems (in addition to being not physically achievable) are not capable of 
showing this behavior because they always implicitly assume an infinite system size 
which effectively discards any interactions from "far away" electrodes. 
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Fig. 4. Polarographic curves for e values of 0, 0.0001, 0.001, 0.01, and 0.1 (listed in order 
from uppermost to lowest curves) with the other physical parameters taken to be 5 = 0, k c = 10, 
and i,. = 10. Notice that for any e > 0, the cell has no problem achieving current densities higher 
than the diffusion limited current (dashed vertical line). All of these curves, with the exception of 
the exact e = curve, were generated by numerically solving Eq. |&U|)- I&61 ) subject to the boundary 
conditions \HT^-I31[l using the method of our companion paper^ 



4.2. Breakdown of the Classical Approximation. For a diffusion-limited 
cell, the classical nonlinear asymptotic analysis just presented leads to an aesthetically 
appealing theory that predicts a limiting current at l = 1. The existence of this 
limiting current fits nicely with our physical intuition that the concentration of cations 
in a solution must always remain nonnegativc. In reality, however, the analysis breaks 
down as the current approaches (and exceeds) its limiting value. 

The breakdown of the classical asymptotics is evident upon examining the expan- 
sions for the bulk field variables as the current is increased toward its diffusion-limited 
value. Calculating a few of the higher order terms in the bulk asymptotic expansion, 



ELECTROCHEMICAL CURRENT- VOLTAGE RELATIONS 19 



we find that 



r«m PM 2l a. 3 2 (2i) 3 , HI 4 (2i) 5 , 6045 6 (2!) 7 8 

(80) - -E(x) = — 7TT- + -e r H — e =• H — e u ttt + 0{e 

V ; c(°) 2 (c(°)) 4 4 (g(0)) 7 4 (c(°)) 



1 £ 2 (2i) 2 3^ (2i) 4 231 c6 (2i) 6 
2 ( g (°)) 2 2 (c(°)) 5 8 (c(0)) c 



(81) c(x) = c(°) + ^ 2 -^3 + ^ 4 7^7^75 + ^^7^778 + 0(e 8 ) 



,2 (2i) 2 6£ 4 J2if_ 777^ (2Q* 
(c(°)) 2 ' (c(°)) 5 4 £ (c(°)) 



(82) p(x)=0 + e — - — 2 +6e // „ N , 5 +-r e TTTTTs + °( e ) 



2 d 2 4> 


e 2 


e 2 


dx 2 


[x + (1 - i)/2i] 2 ' 


X 2 



Since & > — * 2x as l — » 1, the higher order terms are clearly more singular than 
the leading order term at the limiting current. Rubinstein and Shtilman make a 
similar observation from a potentiostatic perspective; they note that the asymptotic 
expansions are not uniform in the cell voltage |24|. 

The inconsistency in the classical approximation was apparently first noticed by 
Levich who observed that the leading-order solution in the bulk predicts an infinite 
charge density when the current density reaches 1, which directly contradicts the 
assumption of bulk charge neutrality 0Sj ■ As l — ► 1, the bulk charge density is given 

by 

(83) p- 

which diverges at the cathode. 

Smyrl and Newman first showed that these paradoxical results are related to 
the breakdown of thermal equilibrium charge profiles near the cathode, leading to a 
significant expansion of the double layer into the bulk solution [53] . They argue that 
the assumption of electroneutrality breaks down when p w c. Since c is proportional 
to x at the limiting current, the bulk approximation fails to be valid for x smaller 
than 0(e 2 ' 3 ), which leads to a boundary layer that is thicker than the usual Debye 
length. From an alternative perspective, the problems begin when p(0)/c(0) ~ 1 
(see Figure |SJ. Using this criterion, we find that the classical asymptotic theory is 
only appropriate when c 3> (2ie) 2 / 3 or, equivalently, i<Cl- (2e) 2 / 3 . Since the cell 
voltage is approximately 4 tanh~ (i) in many situations, this regime also corresponds 
to v — O (| lne|). This shows that in thin films, where e is not so small, it is easy to 
exceed the classical limiting current and achieve rather different charge profiles |22j . 

5. Conclusion. In summary, we have revisited the classical PNP equations, 
analyzing for the first time the effect of physically realistic boundary conditions for 
thin-film galvanic cells and other micro-electrochemical systems. In particular, we 
focus on the effect of Stern-layer capacitance and Faradaic reactions with Butler- 
Volmcr kinetics. Such boundary conditions contain new physics, such as the possibility 
of a reaction-limited current due to the slow injection of ions at the anode. We 
also find that the Stern layer generally allows the cell to exceed limiting currents 
by carrying diverging portions of the cell voltage, which would otherwise end up 
in the diffuse part of the double layer. We have provided analytical formulae for 
current-voltage relations that should provide useful in characterizing the differential 
resistance of thin films, such as those used in on-chip micro-batteries. Here, we have 
focused on the classical nonlinear regime in which thin double layers remain in thermal 
equilibrium; the more exotic, non-equilibrium regime, which arises at and above the 
classical limiting current, is analyzed in the companion paper |22| . 
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Fig. 5. Numerical solutions for the dimensionless cation concentration cj r (x) and (full) 
charge density 2p(x) at the diffusion-limited current (j = l.Oj with physical parameters k c = 10, 
i,. = 10, a c = ct a = 0.5, 6 = 0.0 and e = 0.0001. In the cathode region, electroneutral- 
ity breaks down as the solution becomes cation rich in order to satisfy the reaction bound- 
ary conditions. Note that when x = 0(e 2 ' 3 ), e_|_ and p are both 0(e 2 ' 3 ). For reference, 
the dashed vertical line shows where x = e ' , and the dashed horizontal line shows where 



y = e 2 / 3 [(2 + 2 2 / 3 )+4/(2 + 2 2 / 3 ) 2 ] ^c+(e 2 / 3 ). 
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Appendix A. Positivity of Ion Concentrations . The positivity of the ion 
concentrations follows directly from the mathematical formulation of the problem. For 
the anion concentration, equations 11' I'll can be integrated exactly using the integrating 
factor e^ to yield 



(84) 



c_(.t) = Ae^l 



which implies that the sign of c_(x) is the same across the entire domain. Since 
the integral constraint (|3 1 1> requires that c_(x) is positive somewhere in the domain, 
C-(x) must be positive everywhere in the domain. 

For the cation concentration, we make use of the reaction boundary conditions. 
Integrating equation (|21|l with the integrating factor er, we obtain 



(85) 



c+(x) =c+(0)e 



<K0)-4>(x) 



4ie" 



(*) r 

Jo 



^(y) 



dy. 



Clearly, the integral term is postive because e^ is positive everywhere. Moreover, the 
reaction boundary condition (|29|l implies that c+(0) > because both i r and k c are 
positive. Thus, we find that the cation concentration is strictly positive. 
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